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ABSTRACT 


This  study  investigates  the  effects  of  different  charge  density  rise  times  upon  the 
magnitude  of  the  Cerenkov  radiation  produced  by  a  semi-infinite  electron  beam.  The 
magnetic  field  pulses  were  generated  for  different  rise  times  and  relationships  between 
the  maximum  magnitudes  and  their  associated  rise  times  were  obtained.  These  were 
compared  with  theoretical  relationships  derived  from  a  power  series  approximation.  The 
generated  results  were  close  to  those  predicted  by  theor>'  for  short  rise  times  at  short 
radial  distances.  For  longer  rise  times,  the  departure  from  theory  was  caused  by  the 
magnitude  of  the  decay  portion  of  the  magnetic  pulse.  This  effect  could  be  mitigated  at 
longer  rise  times  by  increasing  the  radial  distance  from  the  electron  beam. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Cerenkov  radiation  is  the  electromagnetic  emission  caused  by  the  passage  through 
a  dielectric  medium  of  a  charged  particle  that  is  traveling  faster  than  the  speed  of  light 
in  that  dielectric  medium. 

Imagine  «Tin  electron  moving  through  a  dielectric  medium  relatively  slowly.  The 
electric  field  of  the  moving  electron  tends  to  repulse  the  electrons  of  nearby  atoms  while 
attracting  the  nuclei,  creating  distortions  in  the  atoms  which  results  in  the  creation  of 
temporary  dipoles.  Thus,  the  medium  around  the  electron  becomes  polarized.  Because 
the  polarization  field  around  the  electron  is  completely  symmetric,  there  is  no  resultant 
field  at  large  distances.  [Ref.  1:  pp.  3-4] 

But  if  the  electron  is  moving  faster  than  the  speed  of  light  in  the  dielectric  medium, 
the  polarization  field  is  no  longer  completely  symmetrical.  While  the  azimuthal  plane 
symmetry  is  preserved,  there  is  a  resultant  dipole  field  along  the  axis  of  travel  which  is 
apparent  at  large  distances.  This  field  is  created  at  each  path  element  along  the  electron's 
track.  Each  field  radiates  a  brief  electromagnetic  pulse  before  the  atoms  realign  them¬ 
selves  and  return  to  their  normal  shape.  [Ref.  I:  p.  4]  This  behavior  can  be  seen  in  Fig¬ 
ure  1. 

The  time  development  of  Cerenkov  radiation  was  introduced  by  Buskirk  and 
Neighbours  for  application  to  electron  beams.  Single  charged  particles  produce  a  radi¬ 
ated  power  that  is  proportional  to  the  frequency.  Electrons  in  an  accelerator  bunch  ra¬ 
diate  coherently,  an  effect  that  more  than  offsets  a  single  particle  increase  in  radiated 
power  with  frequency.  [Ref.  2:  p.  3750] 

This  time  structure  effect  is  not  obser\'able  with  present  technology  in  S  or  L  band 
linear  accelerators  due  to  their  relatively  high  fundamental  frequency.  Induction  accel¬ 
erators  should  produce  observable  results  in  air  for  energies  greater  than  25  MeV  due 
to  their  longer  electron  bunch  structure.  [Ref.  2:  p.  3753] 

Lyman  conducted  a  preliminary  study  investigating  the  magnetic  field  radiated  from 
a  passing  charge  bunch  traveling  over  a  finite  path.  She  showed  that  the  magnetic  field 
magnitude  versus  time  plots  depended  not  only  on  the  observer's  position  but  on  the 
specific  time  conditions  of  each  case.  [Ref  3:  p.  3] 
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Figure  1.  A  Charged  Particle  Moving  Through  a  Dielectric  Medium:  (a)  slow 
speed  (b)  fast  speed  [Ref.  1:  p.  4J 


Price  expanded  on  Lyman's  study  by  investigating  how  the  magnitude  and  shape  of 
the  radiated  field  depended  on  the  observer  location  relative  to  the  time  boundaries  de¬ 
fined  by  relationships  between  the  arrival  times  of  dificrent  parts  of  a  pulse.  [Ref.  4:  p. 
3] 

B.  OBJECTIVES 

One  of  the  key  assumptions  in  the  Cerenkov  radiation  time  development  is  that  the 
current  pulse  has  a  ramp-front,  i.e.  that  the  electron  pulse  density  increases  from  zero 
to  a  maximum  constant  value  in  a  linear  fashion.  The  end  of  the  pulse  is  assumed  to  be 
a  mirror  image  of  the  ramp-front  but  the  cfTccts  of  the  pulse  end  have  not  been  consid¬ 
ered  in  these  studies. 

The  primary  objective  of  this  work  is  to  investigate  the  effects  of  varying  the  rise 
time  of  the  electron  density,  i.e.  varying  the  slope  of  the  ramp-front,  upon  the  magnitude 
of  the  Cerenkov  pulse. 

Secondary  considerations  are  to  compaic  the  shape  and  magnitude  of  Cerenkov 
pulses  along  lines  parallel  to  the  axis  of  electron  beam  travel,  evaluate  the  maximum 


Cerenkov  magnitude  as  a  function  of  radial  distance  from  the  axis  of  electron  beam 
travel,  and  plot  the  magnetic  field  as  a  function  of  time  and  radial  distance. 
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II.  THEORY 


A.  RADIATION  FIELDS 

To  calculate  the  Cerenkov  magnetic  and  electric  fields,  one  must  first  determine  the 
potentials  from  the  moving  charge  distribution.  The  relationship  between  the  fields  and 
potentials  are  given  by  Equations  2.1  and  2.2  [Ref.  2:  p.  3750],  which  are 

B  =  VxA,  (2.1) 

E  =  (2.2) 

Co  dl  ^  ^ 

The  speed  of  light  in  free  space  is  c„  while  the  speed  of  light  in  the  medium  is  c.  The 
electron  pulse  is  assumed  to  move  with  velocity  u  in  the  positive  z  direction.  The  charge 
density  is  assumed  to  be  concentrated  along  the  z  axis  while  the  charge  is  assumed  to 
move  with  no  change  in  shape  [Ref.  2:  p.3750]  such  that 

pM^p(z,s)5(m-),  (2.3) 

PM  =*  Po(^  -  (2.4) 

The  charge  density  p,  has  units  of  charge  per  volume  while  p  and  p„  have  units  of 
charge  per  lenght. 

Using  the  assumption  of  Equation  2.3  and  applying  it  to  the  usual  retarded  solutions 
of  the  wave  equations,  one  obtains  the  potentials  [Ref  2:  p.  3750] 


O(r,0 


i 

e 


■jp{r'/)dz', 


(2.5) 


A(r,/)  = 


p[r',t')dz'. 


(2.6) 


Note  that  c  is  the  permittivity  of  the  medium.  Also,  /'  is  the  retarded  time  and 
R  =  r  -  r'  such  that 
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(2.12) 


The  first  integrals  in  Equations  2. 1 1  and  2.12  have  a  fallofT  of  R-^  at  large  distances 
for  radiation  and  will  be  ignored.  Implementing  the  u  dependence  on  jr  and  j;  (Equation 
2.10  where  culdy=ylR  and  duldx  =  xlR)  into  these  equations  gives  Equations  2.13  and 
2.14  [Ref.  2:  p.  3751J  in  the  form 
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(2.13) 


Po'{u)dz' , 
R 


(2.14) 


Note  that  the  new  term  p/(«)  in  Equations  2.13  and  2.14  is  the  derivative  of  p,(u) 
with  respect  to  the  function  u,  which  was  defined  earlier.  Equations  2.14  and  2.15  can 
be  combined  into  one  equation.  Also  a  new  function  s  can  be  defined  as  5  =  (.x^+^)i 
[Ref  2:  p.  3751]  such  that 


(2.15) 


B.  TIME  DEVELOPMENT 

To  evaluate  Equation  2.15,  one  must  consider  the  dependence  of  u  upon  z‘  as  given 
in  Equation  2.10.  If  Equation  2.10  is  plotted  in  the  m-2'  plane,  the  first  two  terms  are  a 
straight  line  with  unit  slope  and  an  intercept  that  changes  with  time.  The  third  term  is 
a  hyperbola  which  opens  in  the  +  u  direction  with  asymptotic  slopes  of  +  a/c.  In  the 
Cerenkov  case  where  v>c,  the  curve  is  as  shown  in  Figure  2.  The  curve  translates 
downward  to  smaller  values  of  u  as  time  increases.  This  is  due  to  the  time  term  in 
Equation  2.10.  [Ref  2:  p.  3751] 

Certain  assumptions  about  the  charge  shape  and  configuration  must  now  be  taken 
into  account.  The  charge  is  assumed  to  move  with  no  change  in  shape.  This  was  reflected 
in  Equation  2.4:  p{z,t)  =  p,(z  -  vt).  Since  p  and  p,  are  also  functions  of  u,  p{u)  =  p^u). 
The  charge  density  is  assumed  to  have  a  ramp-front  profile,  i.e.  the  density  increases 
from  zero  to  a  maximum  constant  value  in  a  linear  fashion.  The  time  that  it  takes  for 
the  charge  density  to  reach  this  constant  plateau  is  called  the  rise  time  [Ref  4:  p.  13]. 

Only  changing  currents  (those  with  a  nonzero  p,.)  will  contribute  to  the  magnetic 
fields  of  Equation  2.15.  If  the  charge  density  has  a  linear  rise  to  a  constant  value,  then 
the  derivative  p/(w)  will  be  a  constant  valued  square  pulse  of  magnitude  p„'  as  shown  in 
Figure  3.  [Ref  2:  p.  3751] 
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Figure  2.  Tlie  Function  u{z') 


At  large  negative  times,  the  1/(2')  curve  t,  (Equation  2.10)  is  completely  above  the 
nonzero  portion  of  the  p/{u)  curve  so  that  the  contribution  from  p/{u)  is  zero  and  thus 
the  magnetic  field  is  zero.  The  u(2')  curve  moves  downward  as  time  increases  until  the 
B  pulse  begins  when  u(z')  is  tangent  (curve  /, )  to  the  upper  portion  of  the  of  the  /?/(u) 
pulse.  The  magnitude  of  the  magnetic  field  increases  as  m(z')  moves  downward  with  in¬ 
creasing  time  until  the  u(z')  curve  becomes  tangent  (curve  /j)  with  the  lower  edge  of  the 
p,'{u)  pulse.  The  nonzero  part  of  Equation  2.15  has  its  largest  e,xtent  at  this  time-from 
u(zj)  to  1/(23).  At  later  times,  Equation  2.15  breaks  into  two  regions  of  the  2'  axis  and  the 
B  pulse  decreases  with  increasing  time  because  the  extent  of  the  integral  in  the  two  re¬ 
gions  decreases  as  a  result  of  the  upward  turn  of  the  h(z')  curve.  [Ref  2:  p.  3751] 

If  the  derivative  p/{u)  =  p„’  =  constant,  then  it  can  be  taken  outside  the  integral  of 
Equation  2.15.  Also  the  values  of  n  =  eje  and  /?  =  o/c,  can  be  inserted  so  that 


(2,16) 
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The  relationship  of  the  charge  density  p  to  the  function  u  (Figure  4)  must  be  deter¬ 
mined  in  order  to  evaluate  p/.  This  relationship  is 


P 


Po 

A« 


M  +  p, 


0* 


(2.17) 


Note  that  Au  =  m,  -  Wj  (see  Figure  1). 

Taking  the  derivative  of  Equation  2.17  with  respect  to  u,  one  obtains  p„'  such  that 


Po 


/  /  \  /  •  u 

Po  (*^)  ~  Pm  • 

Equation  2.18  is  now  substituted  into  Equation  2.16  to  give 


(2.18) 


B  =  - 


2  Po 
All 


(2.19) 
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Figure  4.  p  as  n  Function  of  u  at  One  End  of  the  Pulse 

Since  *  s’  +  (2  -  z'y  and  Au  —  ua  where  a  is  the  rise  time,  Equation  2.19  can  be 

integrated  over  the  limits  of  2,  to  Introducing  the  function  w  =  z  -  2'  and  integrating 

gives  the  result  , 

» 

®  =  -«/>'  W  ^ -T-  )}■  (2-20) 

Note  that  iv,  =  z  —  z/  and  Wj  =  z  —  z/.  Thus,  the  problem  is  to  find  the  values  of  vv,  and 
w,.  [Ref  3:  p.  14] 

C.  LIMITS  OF  INTEGRATION 

To  find  the  limits  of  integration  z/  and  z/  ,  one  must  consider  two  cases.  The  first 
case  is  where  the  u{z')  curve  minimum  is  in  the  non-zero  portion  of  the  p/(«)  function, 
i.e.  between  curves  r,  and  /j  of  Figure  2.  'I  he  second  case  is  where  the  u{z')  curve  mini¬ 
mum  has  passed  through  the  non-zero  portion  of  the  p,'{u)  function,  i.e.  lower  than  the 
curve  Tj  of  Figure  2. 
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The  first  case  produces  a  single  integral  where  the  limits  of  integration  are  found  by 
solving  Equation  2.10  for  2'.  The  function  11(2')  =  w,  where  m,  is  the  value  of  the  u  curve 
minimum  at  the  upper  edge  of  the  pj{u)  function.  The  result  is 


{^1  -  fz)  ±  VP{(2  -  +  sHl  - 


(2.21) 


Note  that  /I,  =  u,  +  of. 

The  second  case  produces  two  integrals  where  the  limits  of  integration  are  also 
found  by  solving  Equation  2.10  for  2'.  Two  of  the  limits  come  from  Equation  2.21  while 
the  other  two  come  from  where  the  function  u{z')  =  Wj.  t/j  is  the  value  of  the  u  curve 
minimum  at  the  lower  edge  of  the  p/(«)  function.  The  result  is  Equation  2.21  and 


{A,  ~  fi]  ±  +  sHl  -  p'^)) 

(1  -  f) 


(2.22) 


Note  that  /Ij  =  +  of. 

Another  way  to  describe  the  limits  is  that  they  are  the  2'  values  of  the  intercepts  of 
the  u(r')  curve  with  ^he  upper  and  lower  edges  of  the  non-zero  part  of  the  p/(u)  function. 
This  is  shown  in  Figure  5. 

For  the  two  integrals  of  the  second  case,  the  ordering  of  the  limits  is  ver}'  important. 
The  lower  and  upper  linaits  of  the  first  integral  are  z^'i  min)  and  ZjX  min)  respectively 
while  the  lower  and  upper  limits  of  the  second  integral  are  Zii  max)  and  Zi'(  max)  re¬ 
spectively. 


D.  RISE  TIME  DEPENDENCE 

To  derive  an  approximation  for  the  maximum  magnetic  field  dependence  upon  rise 
time,  Equation  2.10  must  be  developed  in  a  power  scries.  Denoting  2/  as  the  value  of 
z'  at  which  Equation  2.10  has  zero  slope,  the  values  of  z/,  u{z,')  ,  and  the  second  deriv¬ 
ative  of  u  with  respect  to  2'  [Ref  2:  p.  3751]  are 


-0  =  2 


•-  1 


(2.23) 


“(V)  =  z-ot  +  s<-Y-l 


(2.24) 
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Figure  5.  The  z'  Limits  of  Integration 


d^u  _  1 
dz'^  ^  0^ 


2D. 


(2.25) 


Z?  is  a  constant  defined  for  conveinence  of  use  Thus  w  can  now  be  expressed  as  a  Taylor 
series  about  the  minimum  sucli  that 

u  =  u{2j)  +  D{2'  -2^'f.  (2.26) 


The  limits  2/  and  Zj'  can  be  written  in  terms  of  the  minimum  value  as  z,'  =  z/  +  Az' 
and  z,'  =  z/  —  Az'  where  Az'  is  the  value  of  z'  —  z/.  The  difference  u(2’)  —  w(z,')  =  a  where 
a  is  the  rise  time  (or  width)  of  the  current  derivative  pulse  pj  (a  can  be  thought  of  as  a 
distance  by  multiplying  it  by  0).  (Ref,  2:  p.  3751]  Therefore 


(2.27) 


Using  this  value,  the  fact  that  Au  =  ua,  and  the  assumption  that  s  and  R  are  slowly 
varjing,  Equation  2.19  can  be  evaluated  to  give 


(2.28) 


Thus  the  maximum  magnetic  field  has  a  dependence  upon  rise  time  to  the  negative 
one-half  power.  Therefore,  as  the  rise  time  increases,  the  maximum  magnetic  field  de¬ 
creases.  [Ref  2:  p.  3752] 

Two  points  should  be  emphasized.  Equation  2.20  gives  the  exact  maximum  magnetic 
field  amplitude  for  evaluated  at  time  ij  (see  Figure  1).  Equation  2.28  gives  an  approx¬ 
imation  for  the  maximum  magnetic  field  amplitude  as  a  function  of  the  rise  time,  giving 
an  opportunity  to  compare  and  contrast  actual  behavior. 
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III.  CALCULATIONS  AND  ANALYSIS 


The  previous  work  had  identified  the  exact  solution  for  the  maximum  magnetic  field 
amplitude  and  an  approximation  for  the  maximum  as  a  function  of  the  rise  time.  Until 
this  time,  no  comparison  had  been  conducted  between  the  two  to  validate  the  approxi¬ 
mation. 

Plots  of  the  magnetic  field  magnitude  versus  time  were  generated  from  evaluating 
Equation  2.20,  using  the  programs  listed  in  Appendix  A.  These  plots  have  200  points 
each.  The  time  increment  between  points  was  adjusted  so  that  the  maximum  magnetic 
field  magnitude  occurred  at  the  fortieth  point  of  the  plot.  The  increment  could  be  so 
adjusted  by  knowing  that  the  maximum  magnitude  would  occur  at  the  end  of  the  rise 
time  a.  The  variables  of  interest  in  the  programs  were  the  rise  time  a,  the  radiil  distance 
s,  and  the  axial  distance  z.  Program  runs  were  made  for  s  =  i.Om,  5.0m,  and  10.0m;  z 
=  0.0m,  1.0m,  and  5.0m;  and  a  =  lO.Ops,  20.0ps,  and  30.0ps. 

Three  general  trends  were  observed  from  these  runs.  As  a  increases  with  s  and  z 
fixed,  the  overall  magnetic  field  magnitude  decreases.  This  is  due  to  the  a->  relation  in 
Equation  2.20.  As  s  increases  with  a  and  z  fixed,  the  overall  magnetic  field  magnitude 
decreases.  This  is  due  to  the  s*’  relation  in  the  arctangent  function  of  Equation  2.20.  And 
finally,  as  z  increases  with  a  and  s  fixed,  the  overall  magnetic  field  magnitude  and  shape 
were  unchanged  except  that  the  pulse  was  delayed  in  time  more  and  more  as  z  increased. 
This  is  because  z  is  merely  a  function  of  time,  as  seen  in  Equation  2.10.  This  time  delay 
behavior  can  be  seen  on  the  horizontal  axes  of  Figure  6,  Figure  7,  and  Figure  8. 

To  further  investigate  plots  of  the  maximum  magnetic  field  magnitude  versus 
the  rise  times  associated  with  each  magnitude  were  generated.  This  was  accomplished 
by  the  programs  listed  in  Appendix  B.  The  magnetic  field  as  a  function  of  time  was  cal¬ 
culated  for  200  points  by  evaluating  Equation  2.20,  just  as  did  the  programs  of  Appendix 
A.  This  was  done  for  a  fixed  s  and  a  fixed  z  (z  was  chosen  to  he  zero  for  all  following 
calculations  since  it  has  no  effect  on  the  magnitude  and  simplifies  the  calculation  proc¬ 
ess).  The  initial  rise  time  was  selected  to  be  lO.Ops.  The  magnetic  field  magnitudes  were 
run  through  a  DO  loop  to  identify  the  maximum.  The  maximum,  along  with  the  corre¬ 
sponding  rise  time,  was  stored  i.-'  an  output  data  file.  The  rise  time  was  then  incremented 
by  5.0ps  and  the  entire  procedure  is  repeated  three  hundred  times.  The  output  data  file 
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Figure  6.  B  versus  r  at  s  =  l.Oin,  z  =  O.Oni,  and  a  =  lO.Ops 
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PULSE  (GAUSS) 


Figure  8.  B  versus  r  at  s  =  l.Oiii,  z  ~  S.Oni,  and  a  =  lO.Ops 
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of  300  points  is  pjotted  by  a  DISSPLA  program  incorporated  into  Appendix  B.  An  ex¬ 
ample  plot  can  be  seen  in  Figure  9. 

The  theoretical  approximation  of  Equation  2.28  predicts  that  the  relationship  be- 
tv’jen  the  maximum  magnetic  field  magnitude  and  the  rise  time  a  should  be  such 
that  oc  0-°^  The  best  uay  to  examine  this  behavior  is  to  apply  a  log-log  scale  to  the 
output  data  file  used  to  create  Figure  9.  If  the  log-log  graph  produces  a  linear  relation¬ 
ship,  the  slope  of  the  graph  will  correspond  to  the  exponent  of  a.  The  log-log  plot 
produces  a  linear  relationship  close  to  that  expected  of  Equation  2.28  for  low  and  me¬ 
dium  rise  times,  but  the  slope  begins  to  decrease  as  the  rise  time  reaches  larger  values. 
This  behavior  can  be  seen  in  Figure  10  and  an  explanation  follows. 

The  exact  solution  for  the  maximum  magnetic  field  magnitude  can  be  calculated  by 
evaluating  Equation  2.20  where  the  limits  z'  (Equation  2.21)  are  for  the  time  (see 
Figure  1 ).  Since  h  =  ti  +  a,  z'  becomes  a  function  of  a.  With  z  arbitrarily  chosen  to  be 
zero,  wjs  becomes  Kajs  where  K  is  some  constant. 

The  next  step  is  to  examine  the  arctangent  function.  For  large  values  of  Kajs,  the 
arctangent  in  Equation  2.20  is  approximately  constant  and  large  changes  in  a  have  little 
eficct  on  the  arctangent  value.  Therefore  the  arctangent  behaves  as  if  there  is  no  a  de¬ 
pendence  and  should  be  proportional  to  a-'  .  But  the  slope  of  Figure  10  is  tending 
towards  zero  for  high  a  .  alucs  instead  of  towards  -1,  the  exact  opposite  of  what  one 
might  think  would  happen.: 

There  is  one  main  reason  why  this  behavior  occurs.;  It  goes  back  to  the  two  ca.scs 
for  the  limits  of  integration  of  Equation  2.20.  The  maximum  magnetic  field  amplitude 
was  calculated  for  the  single  integral  case  where  the  limits  of  integration  were  evaluated 
at  time  (sec  Figure  2).  The  magnetic  field  amplitude  of  the  "tail"  of  the  pulse  was  cal¬ 
culated  for  the  double  integral  case  where  the  limits  of  integration  were  evaluated  for 
times  r>  (see  Figure  5). 

Since  these  calculations  arc  for  the  radiative  fields,  the  magnitude  of  the  magnetic 
field  does  not  decay  to  zero  but  instead  eventually  decays  to  some  constant  value.  As  the 
rise  time  increases  for  a  fixed  observation  point,  the  magnitude  of  the  "tail"  increases 
relative  to  the  maximum  magnitude  of  the  pulse.  Another  way  to  look  at  the  process  is 
that  large  rise  time  pulses  decay  slowly  compared  to  short  rise  time  pulses.  For  extremely 
large  rise  times,  the  maximum  amplitude  is  virtually  indistinguishable  from  the  "tail". 
Note  that  the  maximum  amplitude  for  a  short  rise  time  is  larger  than  the  maximum  rise 
time  for  a  long  rise  time  (from  theory.  oc  Therefore  as  the  rise  time  increases, 
the  maximum  magnetic  field  amplitude  cannot  continue  to  decrease  linearly  because  the 
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maximum  amplitude  is  approaching  a  constant  value  as  the  rise  time  increases.  This 
behavior  is  shown  in  Figure  11  and  Figure  12. 

But  as  s  increases,  it  decreases  the  value  of  Kajs  and  the  arctangent  leaves  the  almost 
constant  value  range  and  tends  towards  the  linear  response  value  region.  As  a  result,  the 
relationship  between  and  a  behaves  linearly  over  a  larger  range  of  a  as  s  increases. 
This  behavior  can  be  seen  in  Figure  13. 

The  data  file  used  in  the  plots  of  Figure  9  and  Figure  10  was  generated  at  a  fixed  s 
value  of  1.0m.  To  further  investigate  this  behavior,  additional  data  files  needed  to  be 
calculated  and  plotted  for  different  values  of  5  to  see  how  much  of  an  effect  the  radial 
distance  might  have  on  the  relationship  between  and  a  .  The  results  are  shown  in 
Table  1. 


Table  1.  VERSUS  a  SLOPE  VALUES 


s  =  1.00m 

slope  =  -0.112 

s  =  5.00m 

slope  =  -0.207 

s  =  10.0m 

slope  =  -0.273 

s  =  15.0m 

slope  =  -0.358 

s  =  20,0m 

slope  =  -0.341 

s  =  25,0m 

slope  =  -0.368 

s  =  30.0in 

slope  =  -0.409 

s  =  40.0m 

slope  =  -0.425 

s  =  50,0m 

slope  =  -0.452 

s  =  75.0m 

slope  =  -0.467 

s  =  100.0m 

slope  =  -0.475 

The  cancelling  of  the  oc  effect  at  high  a  values  by  the  magnitude  of  the 
m'.gnetic  pulse  decay  tail  is  shown  in  Figure  14.  The  effect  of  reducing  the  argument  of 
the  arctangent  function  by  increasing  the  radial  distance  s  and  thus  returning  the 
arctangent  function  more  to  its  linear  response  range  is  shown  in  Figure  15. 

The  relationship  between  the  maximum  magnetic  fip'  '  magnitude  and  the  radial 
distance  s  was  graphed.  The  magnetic  field  was  evaluated  using  Equation  2.20  and  was 
plotted  as  a  function  of  time  for  200  points  at  a  constant  s  and  constant  a.  The  maxi¬ 
mum  magnetic  pulse  magnitude  is  determined  and  stored  with  its  corresponding  s  value. 
The  value  of  5  is  incremented  and  the  entire  procedure  repeated  until  a  300  point  data 
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CERENKOV  PULSE  (GRUSS) 


RISETIME 


Figure  13.  Log  versus  Log  a  at  s  =  3U.0m 
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file  has  been  constructed.  This  data  file  is  then  plotted  using  DISSPLA.  The  programs 
for  these  actions  are  located  in  Appendix  C.  An  example  plot  can  be  seen  in  Figure  16. 

The  last  step  is  to  calculate  and  plot  the  three-dimensional  relation  between  time, 
the  magnetic  field,  and  radial  distance.  The  magnetic  field  is  calculated  versus  time  for 
50  points  using  Equation  2.20.  The  magnitude  and  corresponding  time  and  radial  posi¬ 
tion  are  stored  in  output  files.  The  radial  distance  is  incremented  and  the  entire  proce¬ 
dure  is  repeated  50  times.  The  result  is  three  output  files  containing  the  coordinates  for 
a  50x50  array.  A  3-D  graphing  routine  reads  these  array  points,  converts  them  to  a 
surface,  and  plots  the  surface.  The  programs  for  this  procedure  are  listed  in  Appendbt 
D. 

Due  to  the  fact  that  DISSPLA  will  only  handle  square  arrays  for  surface  plotting, 
the  number  of  points  along  the  B~  i  axis  had  to  be  limited  to  50  to  hold  computer 
memory  requirements  within  acceptable  limits.  This  somewhat  obscures  the  detail  of  the 
3'D  plot.  Therefore  a  limited  2-D  projection  of  B  —  r  curves  for  a  few  5  values  was  cre¬ 
ated  in  Figure  17  to  help  understand  the  workings  of  Figure  18. 

An  example  plot  is  displayed  in  Figure  IS.  If  the  graph  projection  of  the  5  —  s  plane 
is  rotated  about  the  t  axis,  a  representation  of  the  Cerenkov  radiation  cone  in  space  is 
produced. 

Several  values  were  used  in  the  programs  listed  in  the  appendiees.  Electron  energy 
was  set  at  30  MeV.  It  must  be  approximately  25  MeV  or  higher  in  order  for  Cerenkov 
radiation  to  occur.  Current  v.'as  assumed  to  be  510  amperes.  Any  other  desired  values 
can  be  used  for  these  constants.  Except  for  the  purposes  of  Figure  7  and  Figure  8,  the 
observer  axial  distance  r  was  set  to  zero  for  simplification  of  calculations.  The  observer 
radial  distance  s  can  never  be  set  to  zero  except  for  initialization  of  incrementation  use 
in  a  DO  loop.  The  function  w,  w'as  set  to  zero  also  for  simplification.  It  has  the  result 
of  causing  the  function  to  simplify  to  the  expression  =  -  va.  where  a  is  the  rise  time. 
Finally,  the  magnetic  field  calculations  W'ere  conducted  in  .MKS  units  and  converted  to 
CGS  units  from  the  relation:  I  Tesla  «  10,000  Gauss. 
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CERENKOV  PULSE  (GRUSS) 
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B  MRX  VS  RRDIRL  DISTANCE 


RRDIRL  DISTRNCE  (METERS) 


Figure  16.  versus  s  at  a  =  lO.Ops 
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MNITUDE  (GHUSS) 


B-PIELD  VS  TIME  VS  RRDIRL  DISTRNCE 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  primary  objcciive  was  to  determine  the  effect  of  varying  the  rise  time  of  the 
electron  beam  charge  density  upon  the  Cerenkov  pulse  magnitude.  The  theoretical  ap¬ 
proximation  of  Equation  2.28  predicts  that  this  behavior  should  be  oc  a-"-*.  Actual 
results  show  that  this  behavior  approaches  the  theoretical  value  at  large  radial  distances 
due  to  the  decrease  of  the  argument  of  the  arctangent  function. 

This  study  was  conducted  for  the  case  of  a  semi-infinite  path  lenght.  Evaluating  this 
behavior  for  the  finite  case  will  have  some  effect  on  the  expected  outcome. 

The  magnetic  field  magnitude  falls  off  as  the  radial  distance  increases.  This  was  ex¬ 
pected  from  the  relationship  in  Equation  2.20.  Since  the  magnetic  field  must  also  travel 
farther  as  s  increases,  the  pulse  is  time-delayed  relative  to  a  position  closer  to  the 
electron  beam.  By  rotating  the  3-D  graph  of  Figure  10  about  the  time  axis,  an  accurate 
spatial  description  of  the  magnetic  field  pulse  is  obtained.  This  figure  is  the  expected 
Cerenkov  radiation  cone,  somewhat  analogous  to  the  wake  created  by  a  boat  moving 
faster  than  the  current  speed  of  the  water. 

The  shape  and  magnitude  of  the  Cerenkov  pulse  is  not  affected  in  any  manner  by 
changes  in  the  observer's  axial  position.  It  is  only  affected  by  changes  in  the  observer's 
radial  position.  This  assumes  that  all  other  possible  variables  (energy,  rise  time,  current, 
etc.)  are  held  constant.  The  tail  of  such  a  pulse  decays  to  some  constant  value  which  is 
not  zero.  This  effect  is  because  the  calculations  are  for  the  radiative  fields.  Calculations 
should  be  conducted  with  the  non-radiative  fields  as  well. 
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APPENDIX  A.  CERENKOV  PULSE  PROGRAM 


c 

c 

c 

c 

c 

c 
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c 
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c 

c 
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c 

c 

c 

c 

c 

c 


c 


c 

c 

c 


PROGRAM  PULSEI 


*  THIS  PROGRAM  IS  LOCATED  ON  THE  NPS  MAIN  FRAME  COMPUTER  AND  IS  * 

*  WRITTEN  USING  WATFOR  77.  AFTER  ENTERING  CERTAIN  PARAMETERS,  THE  * 

*  SHAPE  OF  THE  CERENKOV  PULSE  IS  GENERATED.  THE  VALUES  OF  THE  * 

*  PULSE  MAGNITUDE  POINTS  FOR  THE  GIVEN  PARAMETERS  ARE  STORED  IN  * 

*  AN  OUTPUT  FILE  FOR  LATER  USE  IN  A  PLOTTING  ROUTINE.  * 


REAL  N,  Ul,  U2,  BETA,  CO,  ROE,  Al,  A2,  BPRME 
REAL  RISETM,  D,  DD,  El,  E2,  Zll,  ZIF,  Z2I,  Z2F 
REAL  Bl,  B2,  Wl,  W2,  YY,  XX,  Tl,  T2,  WPI,  WPF 
REAL  ZPMIN,  SSQ,  BPSQ,  TINC,  lAMP,  ENER,  ENERJ 
REAL  S,  Z,  C,  MO,  V,  NEG,  TPRME,  B 
INTEGER  I 

DIMENSION  TPRME(200),  B(200) 

DATA  TPRME/200*0.0/,  B/200*0. 0/ 

OPEN  CUNIT=4,FILE='TERMINAL') 

OPEN  (UNIT=7,FILE='PULPLT  DATA  A') 


■*Vf**VfVf»WfAyrA5W 


*  N  IS  THE  INDEX  OF  REFRACTION.  lAMP  IS  THE  ACCELERATOR  CURRENT  * 

*  IN  AMPS.  ENER  IS  THE  ELECTRON  ENERGY  IN  MEV.  RISETM  IS  THE  * 

*  PULSE  RISETIME  IN  SECONDS.  S  IS  THE  OBSERVER  RADIAL  DISTANCE  IN  * 

*  METERS.  Z  IS  THE  OBSERVER  AXIAL  DISTANCE  IN  METERS.  CO  IS  THE  * 

*  VACUUM  SPEED  OF  LIGHT  IN  METERS  PER  SECONDS.  MO  IS  THE  ELECTRON  * 

*  MASS  IN  KILOGRAMS.  ROE  IS  THE  LINEAR  CHARGE  DENSITY  IN  COULOMBS  * 

*  PER  METER.  * 

yfVf**iv***Vf*5'r**yf***ivyf*yf****yfyf*Vfyf***yf**vr***Vf*VfVf*yf*yfAVf*yf,M)f****yf******yf* 

« 


N  =  1.000293 
IA.MP  =  510.  0 
ENER  =  30.  0 
RISETM  =  10. OE-12 
S  =  1.  0 
Z  =  5.0 

CO  =  2. 99792458E08 
MO  =  9. 109534E-31 
C  =  CO/N 

ENERJ  =  ENER*1. 6021892E-13 

V  =  CO*SQRT(l. 0-((M0**2. 0)*(CO**4. 0)/(ENERJ**2. 0))) 

IF  (V. LE.C)  GO  TO  60 

BETA  =  V/CO 

ROE  =  lAMP/V 

BPRME  =  N*BETA 


ivyrivyryfivynVivydvycyfAiVjvycyriVivynVivyfyriVAiViVivyfiViVj'ciVAftynvynvyrytivytyrAyfiVivyovyrytyrAyfynv******** 

*  SETTING  Ul  ARBITRARILY  TO  ZERO  FIXES  U2.  THIS  IS  TO  SIMPLIFY  * 
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o  o  o  o  o 


c 

C 

C 


*  CALCULATIONS. 


* 


C 

c 

c 

c 

c 

c 


U1  =  0.  0 
U2  =  -V*RISETM 


*  WITH  200  POINTS,  TINC  WAS  SELECTED  SO  AS  FOR  THE  PULSE  MAXIMUM  * 

*  TO  OCCUR  AT  APPROXIMATELY  20%  OF  THE  GRAPH  LENGHT.  * 


TINC  =  RISETM/40.  0 
SSQ  =  S**2. 0 
BPSQ  =  BPRME**2.  0 
D  =  Z*BPSQ 
DD  =  1  -  BPSQ 
NEG  =  SSQ*DD 


*  ZPMIN  IS  THE  Z  VALUE  OF  THE  MINIMUM  OF  THE  U  CURVE.  * 

»ViV******iV*******A*iViV*'>V'**ilr***Vf****')V*'>Vi*nV>V>V*>VA*****iViV********ilr*Vf**'>V*ft* 


C 

C 

C 

C 

C 

c 


c 


c 

c 

c 

c 

C 

C 


C 

C 

C 

C 

C 

c 


ZPMIN  =  Z  -  (S/(SQRT(BPSQ  -  1.0))) 

***A*jV***5V**j'f*jV*>V*iVA*A******************A*iV********iV>V*iV5V******A**** 

*  T1  AND  T2  ARE  THE  TIMES  THAT  THE  U  CURVE  MINIMUM  IS  TANGENT  TO  * 

*  THE  U1  AND  U2  LEVELS  RESPECTIVELY.  * 

AAAVfAAAiVAAAAyfAVcAAVfAAAAAVfAAAAVfAAAVfAAAAAAAAAAAAAyfAAAAAAAAAAAAAAAAVfAAA 

T1  =  (ZPMIN  -  Ul)/V  +  (SQRTCSSQ  +  (Z  -  ZPMIN)**2.  0))/C 
T2  =  T1  +  RISETM 

DO  40  I  =  1,200 
IF  (I.EQ. 1)  THEN 
TPRME(I)  =  T1 
B(I)  =0.0 
ELSE 

TPRME(I)  =  T1  +  (REAL(I)*TINC) 

END  IF 


AAAAAVfAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAVrAAAAAAVrAAAAA 

*  THIS  IF  STATEMENT  HANDLES  THE  CHANGEOVER  IN  LIMITS  OF  INTEGRAT-  * 

*  ION  FROM  ONE  INTERVAL  TO  TWO  INTERVALS.  * 

AAAAAVfAAAAAVfycyfAAAAAVfVfVfAVfAAAAAAAAAAAAAAAVrAVcAAyfAAAAAAyfAAAAAAAAAAAAAVfA 


IF  (TPRME(I).GT.T2)  GO  TO  10 
A1  =  U1  +  V*TPRME(I) 

El  =  (Z  -  Al)**2. 0 


*  THIS  IF  STATEMENT  HANDLES  THE  CASE  OF  HAVING  A  NEGATIVE  SQUARE  * 

*  ROOT  IN  REAL  SPACE  WHEN  DETERMINING  THE  LIMITS  OF  INTEGRATION.  * 

AAyf.VyfAAAAAyfAyfsVAyfytAyfAAyfAAAAAycAAAAyfAAAAAAAAAyrAyfAAAyfyfAAAAAAAAAAyfyfAAAAA 


IF  (El.LE.ABS(NEG))  THEN 
B(I)  =  0.0 
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GO  TO  30 
ELSE 

Zll  =  ((Al-D)  +  SQRT(BPSQ*(E1+NEG)))/DD 
ZIF  *  ((Al-D)  -  SQRT(BPSQ*(E1+NEG)))/DD 
ENDIF 
C 

WPI  =  Z  -  Zll 
WPF  =  Z  -  ZIF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  ATAN(Wl) 

XX  =  ATAN(W2) 

B(I)  =  -R0E*N*(BETA**2.0)*(YY-XX)’<'1.  0E04/U2 
GO  TO  20 
C 

10  CONTINUE 

A1  =U1  +  V*TPRME(I) 

A2  =  U2  +  V*TPRME(I) 

El  =  (Z  -  Al)*->2.  0 
E2  =  (Z  -  A2)**2.0 
IF  (El.LE.ABS(NEG))  THEN 
B(I)  =  0.0 
GO  TO  30 
ELSE 

Zll  =  ((Al-D)  +  SQRT(BPSQ*(E1+NEG)))/DD 
ZIF  =  ((Al-D)  -  SQRT(BPSQ*(E1+NEG)))/DD 
ENDIF 
C 

IF  (E2.LE. ABS(NEG))  THEN 
B(I)  =  0.0 
GO  TO  30 
ELSE 

Z2I  =  ((A2-D)  4  SQRT(BPSQ*(E24NEG)))/DD 
Z2F  =  ((A2-D)  -  SQRT(BPSQ*(E24NEG)))/DD 
ENDIF 
C 

WPI  =  Z  -  Zll 
WPF  =  Z  -  Z2I 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  ATAN(Wl) 

XX  =  ATAN(W2) 

B1  =  -ROE*N*(BETA**2.  0)*(YY  -  XX)*1. 0E04/U2 
C 

WPI  =  Z  -  Z2F 
WPF  =  Z  -  ZIF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  ATAN(Wl) 

XX  =  ATAN(W2) 

B2  =  -ROE*N*(BETA**2.  0)*(YY  -  XX)*1. 0E04/U2 
B(I)  =  B1  1  B2 
C 

20  CONTINUE 
C 
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30  CONTINUE 

C 

40  CONTINUE 
C 

DO  50  I  =  1,200 

WRITE(7,1000)  TPRME(I),B(I) 

1000  F0RMAT(G15.  7,G15. 7) 

50  CONTINUE 
C 

GO  TO  70 
60  WRITE(4,1010) 

1010  F0RMATC/,'V  IS  LESS  THAN  C,  THEREFORE  THERE  IS  NO  CERENKOV  RADIATI 
$0N  FOR  THIS  ELECTRON’ ,4X,'  ENERGY') 

70  CONTINUE 
C 

STOP 

END 


PROGRAM  PULPLT 

REAL  X(300),Y(300),XMIN,XMAX,yMIN,YMAX 
INTEGER  N 

CHARACTER*20  TITLE$ 

TITLE$  =  'B-FIELD  VS  TIME$' 

N  =  200 

READ(7,8000)  ((X( I) ,Y(I)) , I=1,N) 

8000  F0RMAT(2G15.7) 

CALL  TEK618 

CALL  SHERPA('P6  ','B’,3) 

CALL  PAGE(6,7,5) 

CALL  AREA2D(4.5,5.5) 

CALL  XNAMEC'TIME  (SECONDS)?’ , 100) 

CALL  YNAMEC CERENKOV  PULSE  (GAUSS)?' ,100) 

CALL  HEADIN(TITLE?, 100, 1.5,1) 

CALL  CROSS 

CALL  RANGE (X,XMIN,XMAX,N) 

CALL  RANGEl(Y,YMIN,YhAX,N) 

CALL  GRAF(XMIN,' SCALE' ,XMAX,YMIN, ' SCALE' ,YMAX) 
CALL  CURVE(X,Y, 200,0) 

CALL  ENDPL(O) 

CALL  DONEPL 

RETURN 
END 


SUBROUTINE  RANGE(Y,YMIN,YMAX,N) 
C 

DIMENSION  Y(N) 

YMIN=1. E20 
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YMAX=-LE20 
DO  10  1=1, N 

IF(Y(I).GT.YMAX)  YMAX=  Y(I) 
IF(Y(I).LT.YMIN)  YMIN=  Y(I) 
10  CONTINUE 
C 

RETURN 

END 


SUBROUTINE  RANGE1(Y,YMIN,YMAX,N) 
C 

DIMENSION  Y(N) 

YMIN=1. E20 
YMAX=“1.E20 
DO  20  1=1, N 

IF(Y(I).GT.YMAX)  YMAX=  Y(I) 
IF(Y(I).LT.YMIN)  YMIN=  Y(I) 

20  CONTINUE 
C 

RETURN 

END 
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APPENDIX  B.  B  MAX  -  RISE  TIME  PROGRAM 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


PROGRAM  CERENKV 


*  THIS  PROGRAM  IS  LOCATED  ON  THE  NPS  MAIN  FRAME  COMPUTER  AND  IS  * 

*  WRITTEN  USING  WATFOR  77.  AFTER  ENTERING  CERTAIN  PARAMETERS,  THE  * 

*  SHAPE  OF  THE  CERENKOV  PULSE  IS  GENERATED.  THE  MAXIMUM  VALUE  OF  * 

*  THE  PULSE  MAGNITUDE  IS  DETERMINED  AND  STORED  WITH  ITS  CORRES-  * 

*  PONDIND  RISE  TIME  VALUE.  THE  RISE  TIME  IS  INCREMENTED  AND  THE  * 

*  ENTIRE  PROCESS  IS  REPEATED.  THESE  VALUES  ARE  THEN  READ  TO  AN  * 

*  OUTPUT  FILE  FOR  USE  IN  A  PLOTTING  ROUTINE. 


REAL  N,  Ul,  U2,  BETA,  CO,  ROE,  Al,  A2,  BPRME,  YMAX 
REAL  RISETM,  D,  DD,  El,  E2,  Zll,  ZIF,  Z2I,  Z2F 
REAL  Bl,  B2,  Wl,  W2,  YY,  XX,  Tl,  T2,  WPI,  WPF 
REAL  ZPMIN,  SSQ,  BPSQ,  TINC,  lAMP,  ENER,  ENERJ 
REAL  S,  Z,  C,  MO,  V,  NEG,  TPRME,  B,  DRISE,  BMAX 
INTEGER  I,  J 

DIMENSION  TPRME(200),  B(200),  DRISE(300),  BMAX(300) 

DATA  TPRME/200*0.0/,  B/200*0.0/,  DRISE/300*0. 0/,  BMAX/300*0. 0/ 

OPEN  CUNI’D=7,FILE='PULPLT  DATA  A’) 

OPEN  CUNIT=4,FILE=* TERMINAL') 

OPEN  CUNIT=7, files' CPLOTl  DATA  A') 

*  N  IS  THE  INDEX  OF  REFRACTION.  lAMP  IS  THE  ACCELERATOR  CURRENT  * 

*  IN  AMPS.  ENER  IS  THE  ELECTRON  ENERGY  IN  MEV.  RISETM  IS  THE  * 

*  PULSE  RISETIME  IN  SECONDS.  S  IS  THE  OBSERVER  RADIAL  DISTANCE  IN  * 

*  METERS.  Z  IS  THE  OBSERVER  AXIAL  DISTANCE  IN  METERS.  CO  IS  THE  * 

*  VACUUM  SPEED  OF  LIGHT  IN  METERS  PER  SECONDS.  MO  IS  THE  ELECTRON  * 

*  MASS  IN  KILOGRAMS.  ROE  IS  THE  LINEAR  CHARGE  DENSITY  IN  COULOMBS  * 

*  PER  METER.  * 

*A*»VVf***>V*Vc****A*5Wc************AA**A*A**iV**AA***i^*iViVVf>V5V**?V****'A'*>V** 


N  =  1. 000293 
lAMP  =  510.0 
ENER  =  30.  0 
RISETM  =  9.9E-12 
S  =  10.0 
Z  =  0.  0 

CO  =  2. 99792458E08 
MO  =  9. 109534E-31 
C  =  CO/N 

ENERJ  =  ENER*1.6021892E-13 

V  =  CO*SQRT( 1. 0-((M0**2. 0)*(CO**4. 0)/(ENERJ**2. 0))) 
IF  (V.  LE.C)  GO  TO  80 
BETA  =  V/CO 
ROE  =  lAMP/V 
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BPRME  =  N*BETA 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


*  SETTING  U1  ARBITRARILY  TO  ZERO  FIXES  U2.  THIS  IS  TO  SIMPLIFY  * 

*  CALCULATIONS.  ^  **  **  v* 


U1  =  0.  0 


*  WITH  200  POINTS,  TINC  WAS  SELECTED  SO  AS  FOR  THE  PULSE  MAXIMUM  * 

*  TO  OCCUR  AT  APPROXIMATELY  20%  OF  THE  GRAPH  LENGHT.  * 


SSQ  =  S**2.  0 
BPSQ  =  BPRME**2.  0 
D  =  Z*BPSQ 
DD  =  1  "  BPSQ 
NEG  =  SSQ*DD 


*  ZPMIN  IS  THE  Z  VALUE  OF  THE  MINIMUM  OF  THE  U  CURVE.  * 

*5V****5V*sV***>V**5V5V**5V5V5V*5V***5V»VVfA»VA*5V5V5V5VjV5V***5V****AA5V*5V*)V*A*5VAilriV**5V»V* 


ZPMIN  =  Z  -  (S/(SQRT(BPSQ  -  1.0))) 

****iV**A>VA*iVA**********5V*****AAA**'iV*itiViV****»VAAiV*»VV«VA*7V**iir***A**5V*** 

*  T1  AND  T2  ARE  THE  TIMES  THAT  THE  U  CURVE  MINIMUM  IS  TANGENT  TO  * 

*  THE  U1  AND  U2  LEVELS  RESPECTIVELY.  SINCE  T2  AND  U2  ARE  FUNCT-  * 

*  IONS  OF  RISETM,  THEY  ARE  INSIDE  THE  DO  LOOP.  * 

7V7'fVf7VVfl'f*VflV*5Wr7'C7V)VlVlV5V7V7V7V**lVl'f7V7VAjV7'flV*)VlVlV*7VA7V*lV7V*lV7VlV**7V7VlV***)V7VlV7HtlVsVlV7V*7V*lV 


T1  =  (ZPMIN  -  Ul)/V  +  (SQRT(SSQ  +  (Z  -  ZPMIN)**2.  0))/C 

DO  60  J  =  1,300 

RISETM  =  RISETM  +  1. OE-13 
TINC  =  RISETM/40.  0 
DKISE(J)  =  RISETM 
U2  =  -V*RISETM 
T2  =  T1  +  RISETM 

DO  40  I  =  1,200 
IF  (I.EQ. 1)  THEN 
TPRME(I)  =  T1 
B(I)  =  0.0 
ELSE 

TPRME(I)  =  T1  +  (REAL(I)*TINC) 

END  IF 


***Vt***iV**i’c*****************A*******Vf***Vf****AAVf**yf*************yc** 


*  THIS  IF  STATEMENT  HANDLES  THE  CHANGEOVER  IN  LIMITS  OF  INTEGRAT-  * 

*  ION  FROM  ONE  INTERVAL  TO  TWO  INTERVALS.  * 


iVjvyrsVs’rjvyfyfiVj'tyfi'fAsvyryrsVjViVycyfyfjvyfiViVAyfyfAjVjVi’fiVivyrjVjvyryuVAAAivyfjViVyfyryoVAiVsvycAyryfAAAyfAjViVA 


IF  (TPRME(I).GT.T2)  GO  TO  10 
A1  =  U1  +  V*TPRME(I) 
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El  =  (Z  -  Al)**2.0 


C 

C 

c 

c 

c 

c 


c 


c 


c 


c 


c 


*  THIS  IF  STATEMENT  HANDLES  THE  CASE  OF  HAVING  A  NEGATIVE  SQUARE 

*  ROOT  IN  REAL  SPACE  WHEN  DETERMINING  THE  LIMITS  OF  INTEGRATION. 


IF  (E1.LE.ABS(NEG)>  THEN 
B(I)  =  0.0 
GO  TO  30 
ELSE 

Zll  =  ((Al-D)  +  SQRT(BPSQ*(E1+NEG)))/DD 
ZIF  =  ((Al-D)  -  SQRT(BPSQ*(E1+NEG)))/DD 
ENDIF 

WPI  =  Z  -  Zll 
WPF  =  Z  -  ZIF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  ATAN(Wl) 

XX  =  ATAN(W2) 

B(I)  =  -R0E*N*(BETA**2.0)*(YY-XX)*1.0E04/U2 
GO  TO  20 

10  CONTINUE 

A1  =  U1  +  V’n’PRME(I) 

A2  =  U2  +  V*TPRME(I) 

El  =  (Z  -  Al)**2.  0 
E2  =  (Z  -  A2)**2.  0 
IF  (El.LE.ABS(NEG))  THEN 
B(I)  =0.0 
GO  TO  30 
ELSE 

Zll  =  ((Al-D)  +  SQRT(BPSQ*(E1+NEG)))/DD 
ZIF  =  ((Al-D)  -  SQRT(BPSQ*(E1+NEG)))/DD 
ENDIF 

IF  (E2.LE.ABS(NEG))  THEN 
B(I)  =  0.  0 
GO  TO  30 
ELSE 

Z2I  =  ((A2-D)  +  SQRT(BPSQ*(E2+NEG)))/DD 
Z2F  =  ((A2-D)  -  SQRT(BPSQ*(E2+NEG)))/DD 
ENDIF 

WPI  =  Z  -  Zll 
WPF  =  Z  -  Z2I 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  ATAN(Wl) 

XX  =  ATAN(W2) 

B1  =  -R0E*N*(BETA**2.0)*(YY  -  XX)*1. 0E04/U2 


WPI  =  Z  -  Z2F 
WPF  =  Z  -  ZIF 
W1  =  WPI/S 
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W2  =  WPF/S 
YY  =  ATAN(Wl) 

XX  =*  ATAN(W2) 

B2  =  -R0E*N*CBETA**2-0)*(YY  -  XX)*1.  0E04/U2 
B(I)  =  B1  +  B2 
C 

20  CONTINUE 

C 

30  CONTINUE 

C 

40  CONTINUE 
C 

YMAX  =  -1.0E20 


»v  THIS  LOOP  DETERMINES  THE  MAXIMUM  VALUE  OF  THE  CERENKOV  PULSE.  * 

*****A********^V**A*******************A*****>V********»V************** 


DO  50  I  =  1,200 

IF  (B(I).GT.YMAX)  YMAX  =  B(I) 

50  CONTINUE 

BMAX(J)  =  YMAX 
60  CONTINUE 

DO  70  I  =  1,300 

WRITE(7,1000)  DRISE(I),BMAX(I) 

100  FORMAT(G15.7,G15.7) 

70  CONTINUE 

GO  TO  90 
80  WRITE(4,1010) 

110  F0RMAT(/,'V  IS  LESS  THAN  C,  THEREFORE  THERE  IS  NO  CERENKOV  RADIATI 
$ON  FOR  THIS  ELECTRON' ,4X,’  ENERGY’) 

90  CONTINUE 

STOP 

END 


PROGRAM  CPLOTl 
PROGRAM  LOGPLT 

REAL  X(300) ,Y(300) ,XMIN,XMAX,YMIN,YMAX 
DOUBLE  PRECISION  X(300) ,Y(300) ,XMIN,XMAX,YMIN,YMAX 
INTEGER  N 

CHARACTER*20  TITLE? 

TITLE?  =  'B  MAX  VS  RISETIME?' 

N  =  300 

READ( 7,8000)  ((X( I) ,Y( I) ) , I=1,N) 

8000  F0RMAT(2G15.7) 

CALL  TEK618 
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CALL  SHERPA(*P7 
C  CALL  SHERPAC'PS 

CALL  PAGE(6,7.5) 

CALL  AREA2D(4.5,5.S) 

CALL  XNAMEC'RISE  TIME  (SEC0NDS)$’ ,100) 

CALL  YNAMEC CERENKOV  PULSE  (GAUSS)$’ ,100) 

CALL  HEADIN(TITLE$, 100, 1.5,1) 

PATT  PPHQQ 

CALL  RANGECX,XMIN,XMAX,N) 

CALL  RANGE 1(Y,YMIN,YMAX,N) 

CALL  GRAFCXMIN,' SCALE' ,XMAX,YMIN, ' SCALE' ,YMAX) 
C  CALL  L0GL0G(XMIN,4.S,YMIN,5.5) 

CALL  CURVE(X,Y, 300,0) 

CALL  ENDPL(O) 

CALL  DONEPL 
C 

RETURN  . 

END 


SUBROUTINE  RANGE(Y,YMIN,YMAX,N) 

DIMENSION  Y(N) 

YMIN=1. E20 
YMAX=*1.E20 
DO  10  1=1, N 

IF(Y(I).GT.YMAX)  YMAX=  Y(I) 
IF(Y(I).LT.YMIN)  YMIN=  Y(I) 
10  CONTINUE 

RETURN 

END 


SUBROUTINE  RANGE1(Y,YMIN,YMAX,N) 
C 

DIMENSION  Y(N) 

YMIN=1.  E20 
YMAX=-1,E20 
DO  20  1=1, N 

IF(Y(I).GT.YMAX)  YMAX=  Yfl) 
IF(y(I).LT.YMIN)  YMIN=  Y(I) 

20  CONTINUE 
C 

RETURN 

END 
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APPENDIX  C.  B  MAX  -  S  PROGRAM 


PROGRAM  BVS 

C  . 

C  *  THIS  PROGRAM  IS  LOCATED  ON  THE  NPS  MAIN  FRAME  COMPUTER  AND  IS  * 
C  *  WRITTEN  USING  WATFOR  77.  AFTER  ENTERING  CERTAIN  PARAMETERS,  THE  * 
C  *  SHAPE  OF  THE  CERENKOV  PULSE  IS  GENERATED.  THE  MAXIMUM  OF  THE  * 

C  *  PULSE  MAGNITUDE  IS  DETERMINED  AND  STORED  WITH  ITS  CORRESPONDING  * 

C  *  S  VALUE.  THE  S  VALUE  IS  THEN  INCREMENTED  AND  THE  ENTIRE  PROCESS  * 

C  *  IS  REPEATED.  THESE  VALUES  ARE  READ  TO  AN  OUTPUT  FILE  FOR  LATER  * 

C  *  USE  IN  A  PLOTTING  ROUTINE.  * 

C 

REAL  N,  Ul,  U2,  BETA,  CO,  ROE,  Al,  A2,  BPRME,  YMAX 
REAL  RISETM,  D,  DD,  El,  E2,  Zll,  ZIF,  Z2I,  Z2F 
REAL  Bl,  B2,  Wl,  W2,  YY,  XX,  Tl,  T2,  WPI,  WPF 
REAL  ZPMIN,  SSQ,  BPSQ,  TINC,  lAMP,  ENER.  ENERJ 
REAL  S.  Z,  C,  MO,  V,  NEG,  TPRME,  B,  DS,  BMAX 
INTEGER  I,  J 

DIMENSION  TPRME(200),  B(200),  DS(300),  BMAX(300) 

DATA  TPRME/200*0. 0/,  8/200*0.0/ ,  DS/300*0.0/,  BMAX/300*0.  0/ 

C 

OPEN  (UNIT=4,FILE=' TERMINAL') 

OPEN  (UNIT=7,FILE='BVSPLT  DATA  A') 

C 

C  **Vfyr*****A**Vf**yf*Ayf****Vf****yr*****Vf*A****iVVf**yfilfiV***iViSr**A***Vfyf**Vfilr** 

C  *  N  IS  THE  INDEX  OF  REFRACTION.  lAMP  IS  THE  ACCELERATOR  CURRENT  * 

C  *  IN  AMPS..  ENER  IS  THE  ELECTRON  ENERGY  IN  MEV.  RISETM  IS  THE  * 

C  *  PULSE  RISETIME  IN  SECONDS.  S  IS  THE  OBSERVER  RADIAL  DISTANCE  IN  * 

C  *  METERS.  Z  IS  THE  OBSERVER  AXIAL  DISTANCE  IN  METERS.  CO  IS  THE  * 

C  *  VACUUM  SPEED  OF  LIGHT  IN  METERS  PER  SECONDS.  MO  IS  THE  ELECTRON  * 

C  *  MASS  IN  KILOGRAMS.  ROE  IS  THE  LINEAR  CHARGE  DENSITY  IN  COULOMBS  * 

C  *  PER  METER.  * 

C 

N  =  1. 000293 
lAMP  =  510.  0 
ENER  =  30.  0 
RISETM  =  lO.OE-12 
S  =  0.  0 
Z  =  0.  0 
C 

CO  =  2. 99792458E08 
MO  =  9. 109534E-31 
C  =  CO/N 

ENERJ  =  ENER*1.  6021892E-13 

V  =  CO*SQRT( 1.  0-( (MO**2.  0)*(CO**4,  0)/(ENERJ**2. 0) ) ) 

IF  (V.LE.C)  GO  TO  80 
BETA  =  V/CO 
ROE  =  lAMP/V 
BPRME  =  N*BETA 
C 


41 


o  u  o  o  o 


c 

C 

C 

C 

C 


C 

C 

C 

C 

C 

C 


C 


*  SETTING  U1  ARBITRARILY  TO  ZERO  FIXES  U2.  THIS  IS  TO  SIMPLIFY  * 

*  CALCULATIONS.  * 


U1  =  0.  0 

U2  =  -V*RISETM 


*  WITH  200  POINTS,  TINC  WAS  SELECTED  SO  AS  FOR  THE  PULSE  MAXIMUM  * 

*  TO  OCCUR  AT  APPROXIMATELY  20%  OF  THE  GRAPH  LENGHT.  * 


TINC  =  RISETM/40.  0 
BPSQ  =  BPRME**2.  0 
D  =  Z*BPSQ 
DD  =  1  -  BPSQ 

DO  60  J  =  1,300 
S  =  S  +  0.  01 
DS(J)  =  S 
SSQ  =  S**2.  0 
NEG  =  SSQ*DD 


*  ZPMIN  IS  THE  Z  VALUE  OF  THE  MINIMUM  OF  THE  U  CURVE.  * 

**iV***AVfA**A****iV*iV*iV*5ViV*****5V***AVf**»V**iV*>V*A>V*Vf*iV5V**A***>V*ilr'fr***'iVVc* 


ZPMIN  *  Z  -  (S/(SQRT(BPSQ  -  1.0))) 

C 

Q  A>VAA>V***A*****iV*A*A**iV>VAA***>V*>V*A*AA5V*AA*AA*A**iVA**AAA**A>V**A»WrA*»V* 

C  *  T1  AND  T2  ARE  THE  TIMES  THAT  THE  U  CURVE  MINIMUM  IS  TANGENT  TO  * 
C  *  THE  U1  AND  U2  LEVELS  RESPECTIVELY.  * 

Q  AVcA-.'sAAAAAAAVfVcAAVjVfAAAAVcAAAiVAAAVcAAAAAAAAAAAAAAAAAAAAAAAAAAAVcAAAAAAAAVr 

C 

T1  =  (ZPMIN  -  Ul)/V  +  (SQRT(SSQ  +  (Z  -  ZPMIN)**2. 0) )/C 
T2  =  T1  +  RISETM 
C 

DO  40  I  =  1,200 
IF  (I.EQ.  1)  THEN 
TPRME(I)  =  T1 
BCD  =  0.0 
ELSE 

TPRME(I)  =  T1  +  (REAL(I)*TINC) 

END  IF 
C 

C  *  THIS  IF  STATEMENT  HANDLES  THE  CHANGEOVER  IN  LIMITS  OF  INTEGRAT-  * 
C  *  ION  FROM  ONE  INTERVAL  TO  TWO  INTERVALS.  * 

Q  A******V(A****Ayf*A****Vf*A*************yf*****Vfy{********y(Vc***A**A*Vf*A* 

C 

IF  (TPRME(I).GT.T2)  GO  TO  10 
A1  =  U1  +  V*TPRME(I) 

El  =  (Z  -  Al)**2. 0 
C 
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C  *  THIS  IF  STATEMENT  HANDLES  THE  CASE  OF  HAVING  A  NEGATIVE  SQUARE 
C  *  ROOT  IN  REAL  SPACE  WHEN  DETERMINING  THE  LIMITS  OF  INTEGRATION. 

C  >V**>V*********>V*>V***>V***»V***A******>V*****Vf**iMrjV*>V*********A>V>V*».vV*iHr* 

c 

IF  (El.LE.ABS(NEG))  THEN 
B(I)  =0.0 
GO  TO  30 
ELSE 

Zll  =  CCAl-D)  +  SQRT(BPSQ*(E1+NEG)))/DD 
ZIF  =  ((Al-D)  -  SQRT(BPSQ*(E1+NEG)))/DD 
ENDIF 
C 

WPI  =  Z  -  Zll 
WPF  =  Z  -  ZIF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  ATAN(Wl) 

XX  =  ATAN(W2) 

B(I)  =  -ROE*N*(BETA**2.0)*(YY-XX)*1.0E04/U2 
GO  TO  20 
C 

10  CONTINUE 

A1  =  U1  +  V*TPRME(I) 

A2  =  U2  +  V*TPRME(I) 

El  =  (Z  -  Al)**2.0 
E2  =  (Z  -  A2)**2.0 
IF  (El.LE.ABS(NEG))  THEN 
B(I)  =  0.0 
GO  TO  30 
ELSE 

Zll  =  ((Al-D)  +  SQRT(BPSQ*(E1+NEG)))/DD 
ZIF  =  ((Al-D)  -  SQRT(BPSQ*(E1+NEG)))/DD 
ENDIF 
C 

IF  (E2.LE.ABS(NEG))  THEN 
B(I)  =  0.0 
GO  TO  30 
ELSE 

Z2I  =  ((A2-D)  +  SQRT(BPSQ*(E2+NEG)))/DD 
Z2F  =  ((A2-D)  -  SQRT(BPSQ*(E2+NEG)))/DD 
ENDIF 
C 

WPI  =  Z  -  Zll 
WPF  =  Z  -  Z2I 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  ATAN(Wl) 

XX  =  ATAN(W2) 

B1  =  -R0E*N*(beTA**2. 0)*(YY  -  XX)*1. 0E04/U2 
C 

WPI  =  Z  -  Z2F 
WPF  =  Z  -  ZIF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  ATAN(Wl) 

XX  =  ATAN(W2) 
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B2  =  -ROE*N*(BETA**2.0)*(YY  -  XX)*1. 0E04/U2 
B(I)  =  B1  +  B2 
C 

20  CONTINUE 

C 

30  CONTINUE 

C 

40  CONTINUE 

C 

YMAX  =  “1. 0E2O 
DO  50  I  *  1,200 

IF  (BCD.  GT.  YMAX)  YMAX  =  B(I) 

50  CONTINUE 

BMAX(J)  =  YMAX 
60  CONTINUE 
C 

DO  70  I  =  1,300 

WRITE(7,1000)  DS(I),BMAX(I) 

1000  F0RMAT(G15. 7,G15. 7) 

70  CONTINUE 
C 

GO  TO  90 
80  WRITEC4,1010) 

1010  FORMAT(/,'V  IS  LESS  THAN  C,  THEREFORE  THERE  IS  NO  CERENKOV  RADIATI 
$0N  FOR  THIS  ELECTRON ',4X,'  ENERGY') 

90  CONTINUE 
C 

STOP 

END 

C 

C 

c 

PROGRAM  BVSPLT 
C 

REAL  X(300) ,Y(300) ,XMIN,XMAX,YMIN,YMAX 
C  DOUBLE  PRECISION  X( 300), Y( 300), XMIN,XMAX,YMIN, YMAX 

INTEGER  N 

CHARACTER*25  TITLED 
C 

TITLE$  =  'B  MAX  VS  RADIAL  DISTANCE$' 

N  =  300 
C 

READ(7,8000)  ((X(I) ,Y( I)) ,I=1,N) 

8000  FORMAT(2G15.7) 

C 

C  CALL  TEK618 

CALL  SHERPAC ' P9  ',’B',3) 

CALL  PAGE(6,7.5) 

CALL  AREA2D(4.5,5.5) 

CALL  XNAMEC 'RADIAL  DISTANCE  CMETERS)$' ,100) 

CALL  YNAMEC CERENKOV  PULSE  (GAUSS)?' ,100) 

CALL  HEADIN(TITLE$, 100, 1.5,1) 

PATT  PPnCQ 

CALL  RANGE (X,XMIN,XMAX,N) 

CALL  RANGE1(Y,YMIN,YMAX,N) 

CALL  GRAFCXMIN,' SCALE’ ,XMAX,YMIN, ' SCALE' ,YMAX) 
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CALL  CURVE(X,Y, 300,0) 
CALL  ENDPL(O) 

CALL  DONEPL 

RETURN 

END 


SUBROUTINE  RANGE (Y,YMIN,YMAX,N) 

DIMENSION  Y(N) 

YMIN=1. E20 
YMAX=-1.E20 
DO  10  1=1, N 

IF(Y(I).GT.YMAX)  YMAX=  Y(I) 
IF(Y(I).LT.YMIN)  YMIN=  Y(I) 
10  CONTINUE 

RETURN 

END 


SUBROUTINE  RANGE 1(Y,YMIN,YMAX,N) 

c 

DIMENSION  Y(N) 

YMIN=1.E20 
YMAX=-1.E20 
DO  20  1=1, N 

IF(YCI).GT.YMAX)  YMAX=  Y(I) 
IF(Y(I).LT.YMIN)  YMIN=  Y(I) 

20  CONTINUE 
C 

RETURN 

END 
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APPENDIX  D.  S  -  T  -  B  PROGRAM 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


PROGRAM  STB 


*  THIS  PROGRAM  IS  LOCATED  ON  THE  NPS  MAIN  FRAME  COMPUTER  AND  IS  * 

*  WRITTEN  USING  WATFOR  77.  AFTER  ENTERING  CERTAIN  PARAMETERS,  THE  * 

*  SHAPE  OF  THE  CERENKOV  PULSE  IS  GENERATED.  THE  S  VALUE  IS  INCRE-  * 

*  MENTED  AND  THE  PROCESS  IS  REPEATED.  THESE  VALUES  ARE  STORED  FOR  * 

*  LATER  USE  IN  A  3-D  PLOTTING  ROUTINE.  * 


REAL  N,  Ul,  U2,  BETA,  CO,  ROE,  Al,  A2,  BPRME 
REAL  RISETM,  D,  DD,  El,  E2,  Zll,  ZIF,  Z2I,  Z2F 
REAL  Bl,  B2,  Wl,  W2,  YY,  XX,  Tl,  T2,  WPI,  WPF 
REAL  ZPMIN,  SSQ,  BPSQ,  TINC,  I AMP,  ENER,  ENERJ 
REAL  S,  Z,  C,  MO,  V,  NEG,  TPRME,  B 
REAL  BF,  SRD,  TIME 
INTEGER  I,  J,  K,  L,  M 
DIMENSION  TPRME(200),  B(200) 

DIMENSION  BF(2500),  SRD(2500),  TIME(2500) 

DATA  TPRME/200*0. 0/,  B/200*0. 0/,  M/2500/ 

DATA  BF/2500*0. 0/,  SRD/2500*0.  0/,  TIME/2500*0.  0/ 

OPEN  (UNIT=4, FILE*’ TERMINAL’) 

OPEN  (UNIT=7,FILE='S  DATA  B’) 

OPEN  (UNIT=8,FILE=’T  DATA  B') 

OPEN  (UNIT=9,FILE='B  DATA  B') 

*  N  IS  THE  INDEX  OF  REFRACTION.  lAMP  IS  THE  ACCELERATOR  CURRENT  * 

*  IN  AMPS.  ENER  IS  THE  ELECTRON  ENERGY  IN  MEV.  RISETM  IS  THE  * 

*  PULSE  RISETIME  IN  SECONDS.  S  IS  THE  OBSERVER  RADIAL  DISTANCE  IN  * 

*  METERS.  Z  IS  THE  OBSERVER  AXIAL  DISTANCE  IN  METERS.  CO  IS  THE  * 

*  VACUUM  SPEED  OF  LIGHT  IN  METERS  PER  SECONDS.  MO  IS  THE  ELECTRON  * 

*  MASS  IN  KILOGRAMS.  ROE  IS  THE  LINEAR  CHARGE  DENSITY  IN  COULOMBS  * 

*  PER  METER.  * 


N  =  1.000293 
lAMP  =  510.  0 
ENER  =  30.  0 
RISETM  =  lO.OE-12 
S  =  0.0 
Z  =  0.0 

CO  =  2. 99792458E08 
MO  =  9. 109534E-31 
C  =  CO/N 

ENERJ  =  ENER*1. 6021892E-13 

V  =  C0*SQRT(1. 0-((MO**2.0)*(CO**4. 0)/(ENERJ**2.0))) 
IF  (V.LE.C)  GO  TO  90 
BETA  =  V/CO 
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C 

C 

C 

C 

C 

C 


C 

C 

C 

C 

C 

C 


G 


ROE  =  lAMP/V 
BPRME  =  N*BETA 


*  SETTING  U1  ARBITRARILY  TO  ZERO  FIXES  U2.  THIS  IS  TO  SIMPLIFY  * 

*  CALCULATIONS.  * 

****5VlV*>' *********************************************************** 


U1  =  0.  0 

U2  =  -V*RISETM 


*  WITH  200  POINTS,  TINC  WAS  SELECTED  SO  AS  FOR  THE  PULSE  MAXIMUM  * 

*  TO  OCCUR  AT  APPROXIMATELY  20%  OF  THE  GRAPH  LENGHT.  * 

***5V****************iV******************A********************5V****** 


TINC  =  UISETM/10.  0 
BPSQ  =  BPRME**2.  0 
D  =  Z*BPSQ 
DD  =  1  -  BPSQ 

DO  50  K  =  1,50 

IF  (K.EQ.  1)  THEN 
J  =  0 
ELSE 

J  =  50*(K-1) 
END  IF 
S  =  S  +  0.  5 
SSQ  =  S**2. 0 
NEG  =  SSQ*DD 


*Vf**Vf***yf*Vf******Vf***Vf***Vf*****************************Vf*yf********* 

*  2PMIN  IS  THE  Z  VALUE  OF  THE  MINIMUM  OF  THE  U  CURVE.  * 


v<******vc******Vf********************vc*****Vf****yc***************Vf**** 


ZPMIN  =  Z  -  (S/(SQRT(BPSQ  -  1.0))) 


**Vf****v«*Vf*****v<****Vf**Vf**y.*********Vf**********vt*Vf***************** 

*  T1  AND  T2  ARE  THE  TIMES  THAT  THE  U  CURVE  MINIMUM  IS  TANGENT  TO  * 

*  THE  U1  AND  U2  LEVELS  RESPECTIVELY.  * 

A******A*******ynV**Vf***yf*****A*****A*yr>V**************************** 


T1  =  (ZPMIN  -  Ul)/V  +  (SQRT(SSQ  +  (Z  -  ZPMIN)**2.  0) )/C 
T2  =  T1  +  RISETM 

DO  40  I  =  1,50 

IF  (I.EQ.  1)  THEN 
TPRME(I)  =  T1 
B(I)  =  0.0 
ELSE 

TPRME(I)  =  T1  +  (REALCI)*TINC) 

END  IF 
L  =  J  +  I 
SRD(L)  =  S 
TIME(L)  =  TPRME(I) 
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c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


c 


c 


c 


c 


*  THIS  IF  STATEMENT  HANDLES  THE  CHANGEOVER  IN  LIMITS  OF  INTEGRAT- 

*  ION  FROM  ONE  INTERVAL  TO  TWO  INTERVALS. 


IF  (TPRME(I).GT.T2)  GO  TO  10 
A1  =  U1  +  V*TPRME(I) 

El  =  (Z  -  Al)**2.0 


*  THIS  IF  STATEMENT  HANDLES  THE  CASE  OF  HAVING  A  NEGATIVE  SQUARE  * 

*  ROOT  IN  REAL  SPACE  WHEN  DETERMINING  THE  LIMITS  OF  INTEGRATION.  * 


IF  (El.LE.ABS(NEG))  THEN 
B(I)  =0.0 
GO  TO  30 
ELSE 

Zll  =  ((Al-D)  +  SQRT(BPSQ*(E1+NEG)))/DD 
ZIF  =  (CAl-D)  •  SQRT(BPSQ*(E1+NEG)))/DD 
ENDIF 

WPI  =  Z  -  Zll 
WPF  =  Z  -  ZIF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  ATAN(Wl) 

XX  =  ATAN(W2) 

BCD  =  -ROE*N*(BETA**2.0)*(YY-XX)*1.0E04/U2 
GO  TO  20 

10  CONTINUE 

A1  =  U1  +  V*TPRME(I) 

A2  =  U2  +  V*TPRME(I) 

El  =  (Z  -  Al)**2.0 
E2  =  (Z  -  A2)**2.0 
IF  (El.LE.ABS(NEG))  THEN 
B(I)  =  0.  0 
GO  TO  30 
ELSE 

Zll  =  ((Al-D)  +  SQRT(BPSQ*(E1+NEG)))/DD 
ZIF  =  ((Al-D)  -  SQRT(BPSQ*(E1+NEG)))/DD 
ENDIF 

IF  (E2.LE.ABS(NEG))  THEN 
B(I)  =  0.0 
GO  TO  30 
ELSE 

Z2I  =  ((A2-D)  +  SQRT(BPSQ*(E2+NEG)))/DD 
Z2F  =  ((A2-D)  -  SQRT(BPSQ*(E2+NEG)))/DD 
ENDIF 

WPI  =  Z  -  Zll 
WPF  =  Z  -  Z2I 
W1  =  WPI/S 
W2  =  WPF/S 
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YY  =  ATAN(Wl) 

XX  =  ATAN(W2) 

B1  =  -R0E*N*(BETA**2.0)*(YY  -  XX)*1,  0E04/U2 

WPI  =  Z  -  Z2F 
WPF  =  Z  -  ZIF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  s  ATAN(Wl) 

XX  *  ATAN(W2) 

B2  =  -R0E*N*(BETA**2.0)*CYY  -  XX)*1.  0E04/U2 
B(I)  =  B1  +  B2 

20  CONTINUE 

30  CONTINUE 

BF(L)  =  BCD 
40  CONTINUE 

50  CONTINUE 

DO  60  I  =  1,M 

WRITE(7,1000)  SRD(I) 

1000  F0RMAT(G15. 7) 

60  CONTINUE 

DO  70  I  =  1,M  - 

WRITE(8,1000)  TIME(I) 

70  CONTINUE 

DO  80  I  =  1,M 

WRITE(9,1000)  BF(I) 

80  CONTINUE 

GO  TO  100 
90  WRITE(4,1010) 

1010  FORMAT(/,'V  IS  LESS  THAN  C,  THEREFORE  THERE  IS  NO  CERENKOV  RADIATI 
$ON  FOR  THIS  ELECTRON' ,4X,'  ENERGY') 

100  CONTINUE 

I 

STOP 

END 


PROGRAM  THREED 
C 

REAL  8(2500), T(2500),B(2500),X3MAX,Y3MAX,Z3MAX,ZMAT( 2500) 
INTEGER  I,N 
CHARACTER*35  TITLE $ 

C 

TITLE$  =  'B-FIELD  VS  TIME  VS  RADIAL  DISTANCE?' 

N  =  2500 
C 

READC7,8000)  (S(I),I=1,N) 
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8000  FORMAT(G15. 7) 

READ(8,8000)  (T(I),I=1,N) 

READ(9,8000)  (B(I),I=1,N) 

C 

CALL  RANGE(S,X3MAX,N) 

CALL  RANGE(T,Y3MAX,N) 

CALL  RANGE(B,Z3MAX,N) 

C 

C  CALL  TEK618 

CALL  SHERPAC'PIO  ’,'B',3) 

CALL  PAGE(7.5,9.0) 

CALL  AREA2D(5.  5,7.0) 

CALL  V0LM3D( 1,1,1) 

CALL  X3NAME( 'DISTANCE  (METERS)$' ,100) 

CALL  y3NAME('TIME  (SECONDS)$' ,100) 

CALL  Z3NAME( 'MAGNITUDE  (GAUSS)$' ,100) 

CALL  HEADINCTITLE$, 100, 1.0,1) 

CALL  VUABSCl.  0,-0.5,0.5) 

CALL  GRAF3D(0.  0,' SCALE' ,X3MAX,0.  0, 'SCALE' ,Y3MAX,0.  0, 'SCALE' ,Z3MAX) 
CALL  BGNMAT(50,50) 

CALL  GETMAT(S,T,B, 2500,0) 

CALL  ENDMAT(ZMAT,0) 

CALL  SURMATCZMAT, 1,50,5,50,0) 

CALL  ENDPL(O) 

•;all  donepl 

c 

STOP 

END 


SUBROUTINE  RANGE (Y,YMAX,N) 

C 

DIMENSION  Y(N) 

YMAX  =  ‘1. 0E20 
DO  20  I  =  1,N 

IF  (Y(I).GT.YMAX)  YMAX  =  Y(I) 
20  CONTINUE 
C 

RETURN 

END 
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